Polymorphism and bistability in adherent cells 
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The optimal shapes attained by contractile cells on adhesive substrates are determined by the 
interplay between intracellular forces and adhesion with the extracellular matrix. We model the 
cell as a contractile film bounded by an elastic cortex and connected to the substrate via elastic 
links. When the adhesion sites are continuously distributed, optimal cell shape is constrained by the 
adhesion geometry, with a spread area sensitively dependent on the substrate stiffness and contractile 
tension. For discrete adhesion sites, equilibrium cell shape is convex at weak contractility, while 
developing local concavities at intermediate values of contractility. Increasing contractility beyond 
a critical value, controlled by mechanical and geometrical properties of adhesion, cell boundary 
undergoes a discontinuous transition to a star-shaped configuration with cusps and protrusions, 
accompanied by a region of bistability and hysteresis. 



INTRODUCTION 



Mechanical force generation during cell-matrix adhe- 
sion is strongly influenced by the ability of cells to ac- 
tively probe the mechanical and geometrical cues in the 
extracellular matrix [J. Matrix stiffness plays a pro- 
found role in regulating a variety of cellular processes, 
from morphogenesis, motility to cell spreading and cy- 
toskeletal activity. Cells adhering to softer substrates 
spread less and prefer to have well rounded morphologies, 
while they are more likely to exhibit branched patterns 
on stiffer substrates with greater spread area [2j Ej . Ex- 
periments on micro-patterned adhesive islands revealed 
that cell fate, proliferation and spreading sensitively de- 
pend on adhesion geometry [4]. However, cellular re- 
sponse to extracellular determinants is strongly linked 
to myosin dependent activity of the cell cytoskeleton [5] . 
While myosin activity can influence force transmission by 
regulating the growth of focal adhesions [6], it can also 
drive changes in cell morphology, as seen by pharmaco- 
logically disrupting the cell cytoskeleton [7] or by inhibit- 
ing myosin-II activity [5] . Traction forces exerted by cells 
on substrates can now be determined accurately using 
traction force microscopy or micropillar assays [21 IIP) , 
but the intimate relation between cell shape and trac- 
tion mechanics requires further quantitative investiga- 
tion. In this article we present a minimal mechano- 
geometric model for isolated adherent cells that addresses 
a fundamental question in cell mechanics and morpho- 
genesis: How intercellular and extracellular forces coop- 
erate to control the geometry of cell shapes? At time 
scales when the cell is fully spread and develops stronger 
focal adhesions, the dominant forces in the cell stem from 
surface tension induced by actomyosin contractility and 
elasticity in the actin rich cortex. These intracellular 
forces act in opposition to receptor-mediated adhesive 
forces in determining optimal cell shapes [TTJ [T2]. Al- 
though chemical pathways can trigger a feedback between 
cell activity and cell-substrate adhesion [13], we instead 
focus on their mechanical cooperativity in regulating cell 



shapes. Tuning stiffness of the matrix and acto-myosin 
contractility, we discuss how cells can be driven through 
a series of morphological transitions - convex, concave, 
cusps and protrusions with associated hysteresis. In ad- 
dition, we provide several analytical results relating geo- 
metrical properties of cells e.g. curvature, spread radius 
to mechanical properties such as substrate stiffness and 
contractile surface tension, that are amenable to experi- 
mental verification and quantitative comparison. 

CONTRACTILE FILM MODEL FOR ADHERENT 
CELLS 

We consider a planar adherent cell subject to internal 
contractile forces. The shape of the cell is parametrized 
by the contour r(s), where s represents arc-length. The 
total mechanical energy of the stationary cell is given by: 

E = g J dA + a <j> ds k 2 + k s <j> ds p\r - r Q \ 2 , (1) 

where a is the effective surface tension in the cell due to 
acto-myosin contractility, k is the local curvature of the 
cell boundary and a the associated bending rigidity. The 
last term in Eq. ([I]) represent the strain energy induced 
by the cell on a substrate of stiffness k s through focal ad- 
hesions localized at the cell edge [2] with density p(s), 
so that the total number of adhesions is Na — § dsp. 
The reference configuration r (s) describes the shape at- 
tained by the cell boundary on a rigid substrate, and is a 
constraint set by the adhesion geometry. The assumption 
of local elastic interaction with the substrate is strictly 
relevant for cells cultured on elastomeric pillars or on con- 
tinuous substrates of depth thinner than cell perimeter. 
From a purely mechanical point of view, the Contractile 
Film Model defined by Eq. is equivalent to a fluid 
film bounded by a planar Elastica |15l I16j and attached 
to a continuous or discrete set of links to a compliant 
substrate. Unlike other models of softly restrained films 
p7H2l], Eq. does not involve any constraint on the 
perimeter of the cell, which is then only softly constrained 
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by the adhesion with the substrate. In the special case 
of constant curvature, Eq. ([T]) reduces to the model in- 
troduced in Ref. [7] [221 [23 , with an effective line tension 
A = an 2 . As we will see below, this assumption holds ex- 
clusively when adhesions are uniformly distributed along 
the cell boundary, thus forcing the curvature itself to be 
uniform, while in the case of discrete adhesions k is al- 
ways non-uniform. 

Continuous adhesions. In this case the periphery 
of the cell forms contact with a single continuous adhe- 
sion site, so that p = 1/ C with C — § ds the perimeter of 
the cell. In presence of a uniform and isotropic substrate, 
we can assume the reference configuration to be a circle 
of radius Rq so that a natural minimizcr of the energy 
|l]) would be a circle or radius R. Minimizing Eq. ([!]) 
with respect to R yields the following cubic equation: 



(k s + ira)R 3 - k s R Q R 2 



na 



, 



(2) 



The equation contains two length scales, Rq and £ CT = 
(a/c) 1 / 3 , and a dimensionless control parameter k s /a 
expressing the relative amount of adhesion and contrac- 
tion. For very soft anchoring k s -C a and Eq. (2j admits 
the solution R — £<j, whereas if the cell is rigidly pinned 
at adhesion sites, k s 3> a and R — > Rq. For interme- 
diate values of k s /a the optimal radius R interpolates 
between £ CT and Rq and is an increasing function of the 
substrate stiffness k 8 , in case £ CT < Rq, or a decreasing 
function if £ CT > Rq. For £ CT = Rq, the lower and upper 
bound coincide, and the solution is R = Rq. In par- 
ticular, the case Rq > £ CT reproduces the experimentally 
observed trend that cell projected area increases with in- 
creasing substrate stiffness before reaching a plateau at 
higher stiffnesses [2 |3] . We fit the solution to Eq. Q 
to the measured data for projected areas of 3T3 fibrob- 
lasts adhering to densely spaced pillar arrays of varying 
stiffnesses [24 , as shown in Fig. la. The fitted value for 
surface tension a — 2.7 nN//jm comes to the same order 
of magnitude as reported for endothelial cells [23} [25] and 
epithelial cells [21]. The fit also provides a value for the 
bending rigidity a — 1.2 x 10~ 19 Nm 2 , of the same or- 
der of magnitude as estimated for actin cortex [27] ■ The 
asymptotic behavior and various limits of the solution 
are well captured by the interpolation formula: 



R ■ 



k s R + 3ira £ 
k s + 3ira 



(3) 



indicating that larger surface tension, hence larger cell 
contractility a leads to lesser spread area, consistent with 
the experimental observation that myosin-II activity re- 
tards the spreading of cells [25] • Standard stability anal- 
ysis of this solution under a small periodic perturbation 
in the cell radius shows that the circular shape is always 
stable for any values of the parameters a, k s and Rq. 

Discrete adhesions. For cells adhering to discrete 
number of adhesion sites, one can show that the circular 
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FIG. 1. Geometric properties of an adherent cell, (a) Rela- 
tive cell radius R/Ro as a function of substrate stiffness k s 
for a 3T3 cell cultured on pillar arrays |24j (black). Cell ra- 
dius is estimated from area reported in [24J as R — ^/area/71- . 
Solid (red) line represents the solution to Eq. (J2J with o = 2.7 
nN/pm and a/Rg = 0.06 nN/^m. (b) The total cell length 
£ as a function of adhesion stiffness. For small stiffnesses the 
cell boundary form a curve of constant width (lower inset) 
and £ = nw, with w the width of the curve. This property 
breaks down for larger stiffnesses when inflection points de- 
velops (upper inset). The curvature (c) and the tangent angle 
(d) as function of arc-length for aR^/a = 10, k a Rg = 50 and 
Na = 3. The circle are obtained from a numerical minimiza- 
tion of a discrete version of the energy |l]) , while the solid lines 
are obtained respectively from our analytical approximation. 



solution for the cell boundary is never stable and there 
is always a non-circular configuration with lower energy. 
For simplicity, we assume that Na adhesion sites are lo- 
cated at the vertices of a regular polygon of circumradius 
Rq, with density p(s) — Yli=o 1 — iL), and L the dis- 
tance between subsequent adhesions. Minimal cell shape 
is given by the solution to Euler-Lagrange equation ob- 
tained by minimizing the total energy in ([l]). It reads 



a (2k" + k 3 ) - a + 2k. 



N A -1 

E 

i=0 



5(s — iL) (r — ro) • n = . 

(4) 

Here prime denotes derivative with respect to arc-length 
s and n = r"/|r"| is the normal vector. Due to the 
iV^-fold symmetry of the adhesion sites, adhesion springs 
stretch by an equal amount A in the direction of the nor- 
mal vector: (r, — roi) ■ ni = A, i = 1,2... Na- As a 
consequence of the localized adhesion forces, the curva- 
ture is non- analytical at the adhesion points. Integrating 
Eq. Q along an infinitesimal neighborhood of a generic 
adhesion point i, one finds: 
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FIG. 2. Schematic of cell shapes, (a) <j < <r c i, cell contour is everywhere convex with constant width, (b) a = a p , cell contour 
is purely concave with cusps at adhesion points and without protrusions, (c) a > a C 2, cusps are connected to the substrate by 
means of a protrusion of length I. 



The local curvature of the segment lying between adhe- 
sion points is on the other hand determined by the equa- 



tion, a (2k 



a = 0, with boundary conditions 



n(iL) = n((i + 1)L) = n a . Without loss of generality 
we consider a segment located in s e [0, L], Although 
an exact analytic solution to the nonlinear equation for 
curvature is available (see Ref. [19] and Supporting In- 
formation), an excellent approximation can be obtained 
by neglecting the cubic nonlinearity (Fig. lc-d). With 
this simplification, Eq. Q admits a simple solution of 
the form: 



k(s) = k + — s(s - L) 
4a 



(6) 



Eqs. ([6| and ([5| immediately allow us to derive a con- 
dition on the cell perimeter: L = 2k s A/cr. Furthermore, 
the latter condition leads to a linear relation between 
traction force T — 2k s A, and cell size : 



T = crL 



(7) 



which is indeed observed in traction measurements on 
large isolated cells [25] . 

To determine the end-point curvature kq, we use the 
turning tangents theorem for a simple closed curve [29], 
which requires J Q ds k = 2ir/N A . This leads to following 
relation between local curvature and segment length, or 
equivalently traction force, at the adhesion sites : 
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Finally, to determine the optimal length of the cell seg- 
ment L, we are going to make use of a remarkable geomet- 
rical property of the curve obtained from the solution of 
Eq. Q with discrete adhesions: the fact of being a curve 
of constant width |29) . The width of a curve is the dis- 
tance between the uppermost and lowermost points on 
the curve (see lower inset of Fig. [T|d). In general, such 
a distance depends on how the curve is oriented. There 
is however a special class of curves, where the width is 
the same regardless of their orientation. The simplest 
example of a curve of constant width is clearly a cir- 
cle, in which case the width coincides with the diameter. 
A fundamental property of curves of constant width is 



given by the Barbier's theorem [25], which asserts that 
the perimeter C of any curve of constant width is equal 
to width w multiplied by ir: C — irw. As illustrated in 
Fig. Id, this is confirmed by numerical simulations for 
low to intermediate values for contractility and stiffness. 
With our setting, the cell width is given by: 



w = (Rq- A)(l + cosir/N A ) + h(L/2) , 



(9) 



where h(s) — Jq ds' sinks') is the height of the curve 
above a straight line connecting subsequent adhesions 
and 



0(a) 



ds' k(s') 



+ kos+—s 2 (2s-3L) (10) 
24a 



the angle formed by the tangent vector with the a:— axis 
of a suitable oriented Cartesian frame (Fig. 2a). For 
small angles h can be approximated as : h(s) sa s(L — 
s) [tt/(N a L)~ (<r/48a) s(L - s)}. Using this together 
with Eq. Q and the Barbier's theorem with £ = N A L 
allow us to obtain a quartic equation for the cell length 
and the traction force, whose approximate solution is 
given by: 



crRn 
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9i = 768(1 
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tt 2 )/ [4ttN a (1 + cosir/N A )} and 



Eq. (Ill supports the experi- 



mental trend that traction force increases monotonically 
with substrate stiffness k s before plateauing to a finite 
value for higher stiffnesses [M] [30J. The plateau value 
increases with increasing contractility (Fig. 3a). Trac- 
tion force grows linearly with increasing contractility for 
CTi?p/a -C 1, before saturating to the value 2k s R at large 
contractility aR^/a 3> 1, as shown in Fig.^. Eq. ( fTTj ) is 
also consistent with experimentally observed trend that 
reducing contractility by increasing the dosage of myosin 
inhibitor Blebbistatin, leads to monotonic drop in trac- 
tion forces [3171. 
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INFLECTIONS, CUSPS AND PROTRUSIONS 

For low to intermediate values of a and k s , cell shape 
is convex and has constant width. Upon increasing a 
above a k s — dependent threshold a c i, however, the cell 
boundary becomes inflected (see Fig. |4]and upper inset of 
Fig. [lj>). Initially a region of negative curvature develops 
in proximity of the mid point between two adhesions, but 
as the surface tension is further increased, the size of this 
region grows until positive curvature is preserved only 
in a small neighborhood of the adhesion points. Due to 
the presence of local concavities, the cell boundary is no 
longer a curve of constant width. Convex and concave 
regions are separated by inflection points, given by the 
solution to k — 0, or explicitly: s 2 — Ls + 4aK /a = 0. In 
order for this equation to have real solutions one needs 
ah 3 > 96iTa/N A . Fig. 3d shows <7 c i as a function of k s . 

Upon increasing a above a further threshold value cr C 2, 
the inflected shape collapses giving rise to the star-shaped 
configurations shown in upper right corner of Fig. 4. 
These purely concave configurations are made by arcs 
whose ends meet in a cusp. The cusp is then connected 
to the substrate by a protrusion consisting of a straight 
segment of length £ that extends until the adhesion point 
rest position, so that A ss (Fig. [2j;) (see Support- 
ing Information). The cell boundary becomes pinned 
at adhesion sites as a result of having to satisfy force- 
balance (Eq. Q) and adhesion-induced boundary con- 
dition (Eq. (J5|), while accommodating large contractile 
tensions at its neighbourhood. This results in sponta- 
neous expansion in the cell perimeter. Unlike the pre- 
vious transition from convex to non-convex shapes, this 
second transition occurs discontinuously and is accompa- 
nied by a region of bist ability in the range a p < a < er c2 . 
This is clearly visible in the hysteresis diagram in Fig. 
3c showing the optimal length obtained by numerically 
minimizing a discrete analog of Eq. ([I]) in a cycle and 
using as initial configuration the output of the previous 
minimization. The onset of bistability is regulated by 
substrate stiffness as shown in Fig. [3]i, with stiffer sub- 
strates promoting transition to cusps at lower a C 2- 

At a = o~ p the protrusion have zero length and the 
corresponding configuration of the cell boundary form 
a special curve from which all the shapes having non- 
zero protrusion length can be constructed by means of a 
similarity transformation. To see this let us set I = at 
a = <T p so that the shape of the cell will be of the kind 
illustrated in Fig. 2b. The approximated expression for 
the curvature is the same given in Eq. but with 

A = and «' unconstrained since the last term in Eq. 
Q vanishes identically. The quantities re and L are 
left to determine. A first condition can be obtained by 
observing that: x(L/2) = Rq sirnv /Na, where x(s) is the 
projection of the curve on the edge of the circumscribed 
polygon (see Fig. 2b). A second condition is given by 
the theorem of turning tangents for a simple closed curve 
with Na cusps: J^dsK = 7r(2 — Na)/Na- In the case 
Na = 3, for instance, the right-hand side is equal to 
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FIG. 3. Substrate stiffness and cell contractility control trac- 
tion forces and the stability of cell shapes. Traction force as 
a function of substrate stiffness (a) and contractility (b) ob- 
tained from a numerical minimization of a discrete analog of 
Eq. Q. Solid curves den ote the approximate traction val- 
ues obtained from Eq. (111. (c) Boundary length C obtained 



by increasing (squares) and then decreasing (triangles) the 
contractility for substrate stiffnesses ksR^/a — 100 (green 
squares, black triangles) and k s Ro/a = 120 (red squares, 
blue triangles). The diagram shows bistability in the range 
o p < a < <7 C 2. (d) The critical contractility o c \ and o C 2 as 
functions of substrate stiffness. 



— 7r/3, corresponding to the fact that the tangent vector 
rotates clockwise by 60° as we move counterclockwise 
along the curve from one cusp to the next. These allow 
us to calculate the reference shape shown in Fig. 2b. To 
determine o~ p we use the condition kq = 0, below which 
the curvature turns everywhere negative. An analytical 
estimate for o~ p can be obtained by approximating x(s) 
for small deviations about a flat segment. This gives us: 



3a7r 



N a 



(12) 



Next, following Ref. QH 
balance equation 2k" + k 3 
the scaling transformation: 



Q], we notice that the force 
a / a = is invariant under 



(s, k, a) -> (As, p 



(13) 



Consequently, the equilibrium shape obtained for a given 
value of a > o~ p are similar to the reference shape with 
a scaling factor A = (dp/cr) 1 / 3 < 1. Accordingly, the 
closed curve is rescaled so that L = XL p and A — X 2 A p 
with L p and A p respectively the length and the area of 
the reference shape. This beautiful geometric property 
immediately translates into the following algorithm to 
construct shapes with protrusion (Fig. [2ji) : 1 ) Given the 
surface tension a > o~ p we calculate the scaling factor A. 
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FIG. 4. Cell Polymorphism. Phase diagram in a-k s plane 
showing optimal configuration obtained by numerical mini- 
mization of the energy for Na = 3. 



2) We rescale the reference curve so that L = XL p . 3) 
Finally, we fill the distance between the adhesion points 
and the cusps with straight segments of length I = Ro(l — 
A) (since Rq is the circumradius of the reference shape 
and XRq that of the rescaled shape). This latter step, 
ultimately allows us to formulate a scaling law for the 
length of protrusions that can be tested in experiments: 



£/Ro = 1 - ( V) 1/3 



(14) 



This transition from a smooth shape to a self-contacting 
shape with cusps is reminiscent of the post-buckling sce- 
nario of an elastic ring subject to a uniform pressure |20j . 
But unlike this phenomenon, where the system under- 
goes a continuous transition from a simple curve to a 
curve with lines of contact, here the transition is discon- 
tinuous along both the loading branch (increasing k s ) and 
the unloading branch (decreasing k s ). The transition has 
moreover a strong topological character since it involves 
a jump in the turning number (l/27r) § ds k. 



the Arp2/3 regulated actin array in the lamellipodium. 
Although our model relies on local mechanotransduction 
through adhesions localized at the cell edge, in reality 
traction stresses penetrate inside the cell up to a char- 
acteristic depth controlled by cellular and substrate stiff- 
ness [26] . Our model is thus applicable to cell sizes much 
larger than traction penetration depth and predicts the 
same trend on the dependence of traction forces on sub- 
strate stiffness as derived using long-range elastic mod- 
els [31] . Although, local mechanosensing at cell periphery 
coupled with global surface tension due to cytoskeletal 
contractility can accurately capture experimental trends 
for cell size and traction forces, the effect of non-local in- 
teractions of the cytoskeleton with the substrate cannot 
be neglected at actin remodeling time scales. 

An important consequence of increasing surface ten- 
sion is the loss of stability of smooth shapes and a dis- 
continuous transition to cusps and protrusion (Fig. [3] and 
4 ). The transition is favored on stiffer substrates (see Fig. 
4 ) and leads to spontaneous expansion in the cell perime- 
ter and relaxation of localized adhesion springs. Such a 
transition could also possibly occur on cellular timescales 
via chemo-mechanical instabilities induced by coupling 
of motor activity with ligand-receptor kinetics at adhe- 
sion sites. Instead here it emerges as a consequence of 
the cell boundary satisfying of energy minimization and 
global geometrical constraint imposed by the theorem of 
turning tangents. To our knowledge no experimental evi- 
dence has yet been put forward of such instabilities. The 
result can be tested in cell traction assays by varying 
motor activity in a cycle. 

Finally, recent experiments on multicellular sys- 
tems [55] demonstrated that cohesive cell colonies behave 
like single cells in their distribution of traction stresses 
and presence of supracellular actin cables localized to the 
colony periphery. The Contractile Film Model can be 
conveniently used to study shapes of strongly coupled 
cell colonies, where colony surface tension stems from 
actomyosin contractility as well as strength of cadherins 
mediating cell-cell adhesions. 



MATERIALS AND METHODS 



DISCUSSION 

The Contractile Film Model is limited to the steady 
state cell shapes and does not account for the dynamic 
feedback loop between cytoskeleton and focal adhesions 
during early stages of cell spreading and adhesion. The 
model provides a quantitative framework to describe how 
polymorphic cell shapes arise by tuning substrate stiff- 
ness, adhesion geometry and cell activity. The presence 
of bending deformations in the cell periphery naturally 
allows for optimal cell shapes with non-uniform bound- 
ary curvatures, a feature not included in previous theo- 
retical works on cell shapes [221 [23] . Bending in the 
cell boundary can also arise due to splay deformations in 



Numerical Simulations. The data shown in Figs. 
[T]and[3^,b,d have been obtained by numerically minimiz- 
ing the following discrete version of the energy ([I]) : 

N-l N Na 

Ei = - ^2(x i yi +1 ~Xi +1 yi)+a^2(s l )K^+k s ^2\r i -ro l \ 2 

i—1 i—1 i—1 

where the first term corresponds to the area of the irreg- 
ular polygon of vertices r.; = (xi, y{), with i = 1 2 . . . N, 
and the third sum represents the energetic contribu- 
tion of the Na adhesion points. Ki is the unsigned 
curvature at the vertex i: = |tj — ti_i|/(si) with 
t.- L = (r*i + i — rj)/|r,-+x — ri\ the tangent vector at i and 
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crete energy ( 15 ) was minimized using a standard con- 
jugate gradient algorithm. Using (fl5|) allows a direct 



comparison between simulations and the analytical re- 
sults presented in the previous sections. However, for 
very large substrate stiffness, the discrete curve develops 
self-intersections and the energy becomes ill-defined. In 
this regime, it is more convenient to approximate the cell 
as a simplicial complex consisting of mesh M of equilat- 
eral triangles. The edges of the triangles can then be 
treated as elastic springs of zero rest-length, so that the 
total energy of the mesh is given by: 



The dis- where v and e represent respectively the vertices and the 
edges of the mesh and S is a spring constant. If the 
triangles in the mesh are equilateral, this yields a dis- 
crete approximation of the interfacial energy a A, with 
the spring stiffness proportional to the surface tension: 
i.e. a fa 4273/(2 -B/E), where B/E is the ratio be- 
tween the number of boundary edges B and the total 
number of edges E of the triangular mesh [2TJ . 
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